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ABSTRACT 



The effects of boundary perturbations on finite-amplitude 
acoustical standing waves in rectangular, rigid-walled 
cavities were investigated using non-linear theory. When a 
high amplitude standing wave of frequency w is generated in 
a cavity, non-linear effects v;ill cause a stimulation of 
certain normal modes whose resonance frequencies are integer 
multiples of o). Previous experimental observations revealed 
that there could be excitation of other normal modes, not 
belonging to the family of the driven mode, v;hich v;as not 
predicted by the non-linear theory. 

The purpose of this research v;as to investigate the 
possibility that deviations from the idealized geometry could 
account for these observations. Of the various mechanisms 
possible, this work investigated the possibility of these 
unpredicted excitations occurring through a non-linear 
mechanism. 

The standing waves that exist in an ideal cavity must be 
corrected when the boundaries are irregular. The non-linear 
interaction between these standing waves and the corrections 
was studied. The ability of this interaction to excite 
standing v/aves other than those predicted in the ideal case 
was verified. A specific example was worked out demonstrating 
an unpredicted excitation, the strength of which was on the 
order of the magnitude of the boundary perturbation parameter. 



4 



c 



TABLE OF CONTENTS 



I. INTRODUCTION 9 

II. BACKGROUND 10 

III. THE NON-LINEAR PROBLEM 13 

A. CAVITY DESCRIPTION 13 

B. CLASSICAL SOLUTION FOR IDEAL BOUNDARIES 14 

C. CLASSICAL SOLUTION FOR REAL BOUNDARIES 15 

D. THE APPLICABLE NON-LINEAR EQUATION 16 

IV. THE PERTURBED BOUNDARY 19 

A. DEVELOPMENT OF A BOUNDARY CONDITION ON - 19 

B. A FOURIER SERIES REPRESENTATION 21 

V. A SPECIFIC EXAMPLE 24 

A. THE COSINUS 01 DALLY PERTURBED BOUl^DARY 24 

B. THE SOLUTION FOR P^ 30 

C. RESULTS 31 

VI. CONCLUSIONS 39 

APPENDIX I 40 

BIBLIOGRAPHY 43 

INITIAL DISTRIBUTION LIST 44 



5 



LIST OF SYMBOLS 



Pq Classical linear solution for pressure 

for ideal boundaries 



p. First order perturbation correction due to 
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I. 



INTRODUCTION 



The purpose of this research was to investigate some 
of the effects of boundary wall perturbations on finite- 
amplitude standing waves in a rectangular cavity. The 
investigation was prompted by an examination of the experi- 
mental results of Coppens and Sanders [1] and the research 
of DeVall [2] , which suggested the existence of the 
excitation of modes other than those belonging to the family 
of the driven mode. 
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II. BACKGROUND 



An acoustical wave traveling in a medium will have its 
form distorted as predicted by the non-linear hydrodynamic 
equations. If the medium is absorptive then only waves of 
relatively high amplitude will develop distortion before 
becoming completely attenuated due to the absorptive process. 

At the United States Naval Postgraduate School, Coppens 
and Sanders [1], DeVall [2], and others [3,4] have dealt 
with the study of finite amplitude waves in rigid-v?alled 
rectangular cavities and tubes. One interesting result of 
these cavity experiments was the appearance of excitations 
of modes which v;ere not family members of the driven mode. 

For example, assume that the cavity is driven in the (0,1,0;1) 
mode, which indicates the excitation is of the form 
cos k^y cos wt , where 



and 




Then a non-familial mode is any which does not have the 
form (0,n,0;n) . It was proposed [5] that the presence 

of these non-familial modes could be attributed to the cavity 
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walls being non-ideal (not perfectly flat) and that perhaps 
the previous small disparity between theoretical predictions 
and experimental observations might be partly explained by 
the presence of these non-familial modes. 

The non-linear wave equation applicable to this problem 
can be expressed [6] as 



2^2 p 

^0 Dl : 



19^ r / P ^2 . ,^0^2, 

2 — 2 2 ^ ^ 

3t2 Cq 



( 2 . 1 ) 



^ 1 ^ 2 „2 r / P N 2 2, 

+ 2 Cq V [( ^) - (— ) ] 




is defined as 



□ l = D" - L ' 

where L represents the loss terms appropriate to each 
frequency component in P . The RIIS of the non-linear 
equation may be interpreted as a forcing function (consisting 
of both classical and non-linearly generated terms) which, 
in first-order perturbation theory, is the classical solution 
for the pressure. The first-order perturbation solution may 
be considered to describe the effects of self-interaction of 
the classical solution, whereas higher-order perturbation 
solutions describe the effects of interaction of the 
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of the non-linear solution with itself. If the cavity is 
driven at frequency w , non-linear terms force the 
existence of all terms containing the form nw . Each such 
term whose frequency is close to the resonance frequency of 
a standing wave of the cavity and whose spatial behavior 
matches that standing wave can strongly excite that standing 
wave to a degree determined by Q , the quality factor of 
the resonance, and by the difference between nw and the 
resonance frequency of the standing wave. In a rectangular 
cavity, it can be seen [1] that if the classical solution 
is the (m ,m ,m ;q) mode, then the non-linear wave equation 
will force the strong excitation of all modes 
(nm„,nm_,nm^;nq) which constitute the family. Any mode of 
frequency nw which is not a member of the family should 
not be strongly excited. The purpose of this research is to 
see if the presence of wall irregular! tes can explain how 
non-family members may be strongly excited, and to present 
an example to support this theory. 
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III. THE NON-LINEAR PROBLEM 



A. CAVITY DESCRIPTION 

Assume a perfectly rigid-walled rectangular cavity has 

one wall perturbed such that the cavity dimensions are 

L [l+6f(y)l, L , L . A top view depicting the XY plane 
X y z 

and showing the boundary irregularity is given in Fig. 1 
below. Also assume the perturbation on the boundary is 
small compared to the cavity dimensions, l6f(y) ] << 1 . 




FIGURE 1 



The cavity is to be excited by a source near one of the 
corners in such a way that the (n,m,5.;q) mode is excited. 



13 



B. CLASSICAL SOLUTION FOR IDEAL BOUNDARIES 



For a rigid-walled rectangular cavity with ideal 
boundaries (6=0) , the classical linear solution for the 
pressure is subject to the following conditions, 

namely 



□ % = 0 (3.1) 

and Vpp'n^ =0 at x=0,L^ (3.2) 

y=0 ,L 

y 

z— O/L^ / 

where n^ is the local normal to the ideal surface. The 
solution for p^ in terms of Mach number is 



= M cos k X cos k y cos k z cos(co t+cf)) , (3.3) 






0^0 



where 



X 



mr 

L 

X 



(3.4) 



miT 

L.. 



g,TT 

L 
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and 




( 3 . 5 ) 



where 03 is the frequency associated with the (n,m,i!.;q) 
Si 



mode. 

C. CLASSICAL SOLUTION FOR REAL BOUNDARIES 

Since the cavity is assumed to have one perturbed wall, 
it is reasonable to expect the solution for real boundaries 
to differ from that for ideal boundaries. The classical 
linear solution for real boundaries, , may then be 

considered a summation of the classical linear solution for 
ideal boundaries plus perturbation correction terms due to 
the irregular boundary: 



As long as the magnitude of the boundary perturbation is kept 
small, terms higher than first-order can be considered 
insignificant, so that 



P 



0 




P 



0 



Pq + 6pj^ (to first order) 



( 3 . 6 ) 



and must be a solution to 




( 3 . 7 ) 
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with boundary conditions 



VP^-n = 0 at x=0, L^[l + 6f(y)] (3.8) 

y=0 , L 

y 

z=0, L , 
z ' 

/s 

where n , the local normal to the real surface, is obtained 
by taking the divergence of the equation for the perturbed 
boundary, i.e., V • {x - [l+6f (y) ] } . Thus, to first order 

in 6 , 



n 



X 






3f (y) 

9y 



y 



f 



and when the above is used in (3.8) the result is 



3P 



0 



9x 



x=L^[l+6f (y) ] 






^^0 9f (y) 
9y 9y 



x=L^[l+6f (y) ] 



= 0 . 
(3.9) 



D. THE APPLICABLE NON-LINEAR EQUATION 

A real wave in a real cavity will of course have non- 
linearities introduced as it propagates through the cavity. 

Let P^ be the correction term to the pressure due to the 
non-linearities of the medium. The total pressure P is tlien 

P = Pq ^1 (3.10) 
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I'Jhen (3.10) is inserted in (2.1) , the LHS of the non-linear 
equation becomes 



r. 2_2 

0 Ql 






0 ^ 2^2 
— + ^0 Ol 



Po^o 



The first term above is known to be zero from (3.7) , so the 
LHS of (2.1) becomes simply 



2^2 
*^0 Ol 






Now if standing waves are to be produced in the cavity, P 



0 



must consist of terms of the form cos k^x cos k^y cos k^z cos oit 
Using the basic definitions 



34.Q 
‘’O 8t 



and 







0 



requires that U,. be of the form sin k x sin k y sin k z sin wt 

u X y z 

When this information is used, the RHS of (2.1) will consist 
of four terms : 



a 2 

(cos k X cos k y cos k z cos tot) 
at'^ ^ ^ ^ 
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2. 



9^ 2 

— TT (sin k X sin k y sin k z sin wt) 

8t ^ ^ ^ 



2 2 
3. V (cos k X cos k y cos k z cos wt) 

X y z 



2 . 2 
4. V (sin k x sin k y sin k z sin wt) 

X y z 



Although the expansion of these will produce many terms, it 
must be noted that only some will produce resonant standing 
waves. All the other terms may therefore be considered to 
give a small contribution when compared to the resonant terms, 
and thus may be neglected. This fact, when used with the 
trigonometric identities 



and 



2 1 

cos a = 2 [1 + cos 2 a] 

2 1 

sin a = 2 tl“COS 2 a] 



reduces the RHS of (2.1) to a single term, and after a good 
deal of manipulation the non-linear equation may be 
approximated by 



^22 

^0 Dl 



Y + 1 y 



^)2 






9t' 



PnC 



0"'0 



(3.11) 



The governing non-linear equation is now in a relatively 
simple form, and Appendix I develops the general form of 
the solution for . 
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IV. THE PERTURBED BOUNDARY 



A. DEVELOPMENT OF A BOUNDARY CONDITION ON 

Recall that is the first-order correction term due 

to one boundary being irregular. Although it is known that 
for a perfectly rigid-walled cavity the pressure must be a 
maximvim at the walls , this fact may not be conveniently 
used with one wall having an irregular form such as 
L [l + 6f(y)] . It would be highly desirable, therefore, 

to develop a boundary condition on p, to be applied at L . 
This may be accomplished as follows; 

Using the incremental form for a Taylor series expansion 
on P» evaluated at the real boundary L [l+6f(y)] produces 

U X 



x=Lx[l+<Sf (y) ] 






8x 



X=L 



x=L 



Vf(y) + 

3x 



g^p 

1 o 2 

^ ^ [L^6f(y)]^ + ... 



X 



x=L. 



Substitution of 



Pq = Pq *^Pl 

on the RHS, and taking the partial derivative with respect to 
X on both sides gives 
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8x 



x=Lx[l+6f (y) ] 



9P( 

9x 



+ 6 



x=L 



9P] 

9x 



X 



x=L. 



X 



9 Pr 



3x 



x=L 



L^6f(y) + 



(4.1) 



X 



Following exactly the same procedure for the y coordinate 
gives 



9y 



y=Lx[l+6f (y) ] 



aP( 

9y 



+ 6 



x=L 



9P] 

w 



X 



x=L 



X 



^ Pq I 
8y^ lx=L 



L 6f (y) + . . . . (4.2) 

X 



X 



Recall that (3.9) required 



9P 



0 



9x 



x=L [l+6f (y) ] 






^^0 9f (y) 
9y 9y 



x=L [l+6f (y) ] 

X 



0 . 



When (4.1) and (4.2) are used in (3.9) the result is 
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8P( 

W 



6 8P- 



x=L 



9x 



a^p, 



X 



x=L^ 8x 



x=L 



L^6f(y) + 



X 



x°ay ^3y 



63Pi 



x=L 



3y 



3^Po 



3y 



X 



x=L 



L^6f (y) + , 



. ]= 0 , 



X 



Collecting like-order terms in 6 then gives 



9Pi 



3x 



9^P, 



x=L 3x 

X 



L f (y) - L ^ 

X ^ X 3v 9v 



x=L 



X 3y 3y 



= 0, 



X 



X=L 



X 



or, as a boundary condition on applied at L , 

_L X 



ap. 

ax 



9^P, 



X=L. 



= L t-f(y) _j 



X 



8x 



x=L. 



af (y) 
ay ay 



X 



]. (4.3) 



x=L 



X 



B. A FOURIER SERIES REPRESENTATION 

It will be desirable for further development to represent 
the RHS of (4.3) as a Fourier series in cosines, so that p^ 
may be expressed as a summation of normal modes. First, 
rewrite (4.3) as 



ax' 

Po^o 



jML^[-f (y) 



x=L 



2 

3 Po 

7 

ax 



x=L 



af (y) 
ay 



ay 



(4.4) 



x=L 



X 



X 
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Let the RHS of (4.4) equal Scos k^z cos cot and express it 
as an infinite Fourier series as follows: 



S(y)cos k^z cos ojt = [-— + 2 cos cos k^z cos wt, 

N=l^ y 



where the are given by 



L 

2 - y 



= jf- ^ S(y) cos ^ZZ dy, 
y 0 y 



N = 0,1,2, 



Then (4.4) becomes 



J-( “"i 1 

3x' 



P c 
0 0 x=L 



Aq 00 

= (- 2 “ + E Aj^ cos -^) cos k^z cos o)t. 
N=1 y 



The solution for is 




P c 
^0 0 



2 A COS cos k z cos wt (4.5) 

N k sxn k L y 

XXX ^ 



where for N = 0 the coefficient is A„/2, and k.' is such that 

U 



k 



I 

X 



/(f," - 

s! y 



Note that k^' is not in the form of a standing wave in the 
cavity required by 
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(3.4) 




Thus the terms involving (4.5) must also be expressed 

as a Fourier series: 

cos k 'x 

X 

k ' sin k ' L 

X XX 

where the are given by 

L u I 

P X cos k 'x 

= j— / cos dx, M = 0,1,2, .. . 

^x 0 k^' sin kj^'L^ ^x 

When the results of (4.6) are inserted into (4,5), is 
expressed as a double infinite series: 



. ^ 



B 



M=1 



M 



cos 



Mttx 



(4.6) 



^ = - E Z Bj^ cos cos cos k z cos cot, (4.7) 

PqCq N M X y ^ 

where again each of the zero coefficients (A , B ) are to be 

0 0 

divided by a factor of two. 

The formalism has now been developed whereby given any 
driving mode and boundary irregularity, p^ - the first-order 
perturbation correction term - may be expressed in terms of 
the spatial components of the normal modes of the ideal 
cavity. 
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V. A SPECIFIC EXAMPLE 



A. THE COS INUSOI DALLY PERTURBED BOUNDARY 

As a specific example utilizing the theory developed 
in the preceeding sections, assume the cavity is driven in 
the (0,1,0;1) mode. Thus 






= h{ cos k y cos wt 
^ y 



(5.1) 



Further assume that the boundary is perturbed cosinusoidally 
such that 



f (y) = 



JL 

L 

y 



(5.2) 



The boundary condition on p^ then becomes 



^ 









x=L 



V . TTV 

W = — TT k Sin ~ sin k„y cos mt 
L y L y' 

y y 



X 



which, when used with the requirement that k 
becomes 



TT 



y L ' 
y 



3 , ^1 









x=L^ 



M ~ k cos mt - M k cos cos wt 

Ly 2 y Ly 2 y Ly 



X 
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By inspection, the solution for must be 



PnC 



0"'0 



L k cos k 'x 

“ IT 2 iT^ sin k 'L 
y X XX 



L k cos k "x - 

, ,, X IT y X 2TTy 

^ L 2 iT^ sin k "L ir~ 

y X X -^x y 



(5,3) 



where k^' satisfies the equation 



X c 



(5.4) 



and k^" satisfies 



( k ^")2 + (k + = (^)2 

y 



(5.5) 



Solution for p^ allows the creation of Pq = p^ + 6p^^ , 



PQ^O 



L k cos k 'x 

M cos k y cos cot - 6M rr ■ ■■ ■ I V ■ cos tot 

y L 2 k ' sin k 'L 

y X X X 



j xi/ ^x IT ^y ^x ^ 2iry 

+ — 2 ir^ k ■*" L • 

y X X y 



(5.6) 



Let the cavity dimensions be designed such that the (0,2,0; 2) 
and the (1,1,0; 2) modes will be degenerate. This requires 
that 
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(2 






and 



(ir> ' - <iT) ' 

X V 



, 2u) 2 

' c ^ 



Solution of the above equations yields 




(5.7) 



for the desired degenerate modes. t'Then (5.7) is inserted 
into (5.6), the result is 






= M cos k y cos cot - 6 W • 5 - 



i ^ cos 

/r 2 *^x' 



COS cot 



+ 6 M 



nr 2 k " 
/3 X 



cos k X 

X 2iTy . 

— ; ^ irr- cos cos COt 

sin k "L L 

XX y 



(5.8) 



This must be squared when substituted into (3.11). 
Recalling that 



P 



0 



Pq + «SPi 



(3.6) 
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then the non-linear term has the form 



^0 2 Pq 2 
( %) = ( + 26 (- 









p c 
^0 0 



r) ( 



'1 



r) + <5 ( 



2, Pi 



•) 



~ 2 > • - ' 2 ' . 

P «c pc 

0 0 0 0 



(5.9) 



2 

The third term on the RHS above is of order 6 and may 
therefore be neglected. The first term becomes 



, Pq X ^ ,,2 2 , 2 ^ 

( = M cos k y cos wt. 

2 y-' 

pc 
^0 0 



Using the trigonometric identity 



2 1 

cos a = •^ [1 + cos 2a] 



produces the following for : 



Pn 2 ..2 

( [1 + cos 2k y] [1 + cos 2wt] 

^0 0 



Since the only significant contribution will come from 
terms near resonance, those terms not involving a product 
of cos 2k^y cos 2cot may be neglected. Thus 






PO 2 . ^2 

2^ 4 



cos 2k y cos 2wt 

y 



(5.10) 



27 



similarly, keeping just those terms which will lead to 
resonance produces 



26 



^0 



Pt j.1 o o 1 k COS k ' X 

1 _ _ Y+1 1 IT y X 

"2" 2" ~ ~ "2 — ^ ^ ( 2 (jo ) — "2 ^ — 



^00 ^00 



- cos k ycos 2 (jjt 
k ' sin k 'L ^ 

X XX 



o o 1 ^ cos k "x 

+ -i- 5 — 6 M (2u)) ^ 5 - 

/J ^x sin k "L 



X X 



cos k^y cos 2 u)t . 

(5.11) 



Inserting the results of (5.9 - 5.11) into (3.11) gives 



2_2 1 Y +1 .,2 2 

^0 Ol 1 ~ 2~ ^^y^ 

PQ^O 



o o T k COS k 'x 

Y+1 ,,2 ^2 1 -IT y X 

- 6 M (2w) rr — ^ 

^ r=r 2. 



- cos k y cos 2 'a)t 
/3 "• k. ! sin k. 'L ^ 



X X 



,1 o o T k COS k "x 

+ 6 (2a,) 2 A f 2 

^ ’'x" *'x"^x 



cos k y cos 2 o)t. 

(5.12) 



The terms cos k^'x and cos k^"x must be expressed as a cosine 
series of the form given by (4.7) . When this is done, and 
only those components of argument ttx/L are retained, then 
all terms driving the ( 1 , 1 , 0 ; 2 ) mode will be isolated, as 
desired. The results are 



k 



I 

X 



= TT/Ly 



(5.13a) 
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cos k X = 0.535 cos = — , 
X L 

X 



(5.13b) 



sin k 'L = 0.971 , 



(5.13c) 



k^" = i /3~ IT 



(5.13d) 



cos k "x = 



-3.677 cos 



TTX 

L ' 
X 



(5.13e) 



and 



sin k "L^ = i23.1 

X X 



(5.13f) 



where use was made of the identities 



sin iA = i sinh A 



and 



sinh A = •j [e - e ] . 

When the results of (5.13) are inserted into (5.12), the 
result is 

^O^Dl ^ 

- 61.83 cos ^ cos ^ cos 2cot, (5.14) 
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where the first term on the RHS is called the finite 
amplitude term, and the second term on the RHS shows the 
effects of the boundary irregularity on the first order 
non-linear perturbation correction. 

B. THE SOLUTION FOR 

Appendix I shows the general solution for a function 
2 2 

operated on by c . Using those results, the solution 

for for this example becomes immediate: for the finite 

amplitude term, designated by subscript f, 

^1 Y+1 ^f 2 2 

( ^ M u cos 2k y cos (2wt+3^.) ,(5.15) 

PqCq f (2u) ^ 



and for the boundary irregularity induced term, designated 
by subscript p. 



( 



QqCq P 



-1.83 







2 Q sin B 

co^ _E cosk xcosk ycos (2wt+B„) • 

(2w)^ X y P 

(5.16) 



When (5.15) and (5.16) are combined, they form P^ , the first 

order non-linear perturbation correction to the pressure. 

P, is a non-linearly generated term which results from 
■^f 

the existence of a 010 standing V7ave in the cavity. P, 

P 

exists because the irregular boundary forces a non-linear 
interaction between the correction term p^ and the 010 mode. 

It should be noted that there exists an additional 
correction term to the total pressure. This is the classical 
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correction to the 020 standing wave which exists in a 

cavity with an irregular wall. This additional correction 

term will be on the same order of magnitude as [7] , 

P 

and must be taken into account when finding the total effect 
of the irregular boundary on the pressure. The following 
discussion of results will consider only the correction term 
generated by this example. 

C. RESULTS 

In order to best illustrate the magnitude of the effect 
of the non-familial mode in this example, recall that from 
(5.15) and (5.16) , 



^1 Y+1 ^f ^f 2 2 

2 ' ~ 2 — 2 — ^ ^ (2wt+3£) 



Po^o 



( 2to) 



^ Q sin 3 9 ? 

- 1.83 5 ^ 9 —^ M u cos k x. cos k y cos ( 2 cot +3 ) 

(2w) X y p 



Dividing both sides of the equation by 



Y+1 
2 4 



allowing the quality factors and Q to be equal, and 

**" ir 

assuming measurements are made at one of the corners so that 



cos k X = cos k y = cos 2k y = 1, 

X y-^ y 
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gives 



P,/p c 
inn 

a±i £ Q 

2 ^ 



sin 3^ cos (2wt+3£) 



1 . 836sin 



3 cos 
P 



(2wt+3p) 



Let the LHS be designated by C. Additionally, introduce 



H 



Q,. Q 

f p 

Qx: + Q 

p 





as a measure of the overlap of the curves obtained by 
plotting the magnitude of C versus a function of frequency 
for two cases : 

1 - 6=0, corresponding to no non-familial mode 

existence, and 

2 - 6 0 , corresponding to a degeneracy between the 

(0,2,0; 2) and (1,1,0; 2) modes. 



These curves. Figs. 2-6, allow a comparison of the magnitude 
of the term with and without the presence of a non- 
familial mode. The function of frequency on the abscissa 
is 



F = cot 3j^, 

where 3 is as defined by (I.l) . Plots are provided for five 
different values of H. Additionally, for the case H = -1 
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FIGURE 
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FIGURE 




35 



FIGURE 
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FIGURE 
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FIGURE 



there is a plot labelled H_ = -1 which shows the result of 
selecting a different corner for measurements, one 180° 
out of phase with the plot labelled = -1. 6 was arbi- 

trarily chosen as 0.1. 

From the graphs it is clear that the presence of a non- 
familial mode does affect the total solution for the pressure 
by an amount proportional to 6, the magnitude of the boundary 
perturbation . 
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VI. CONCLUSIONS 



Non-linear theory has been applied to standing waves 
in a rigid-walled rectangular cavity with a perturbed 
boundary in an attempt to find one possible mechanism for 
the excitation of resonant modes other than those belonging 
to the family of the driven mode. It was observed that such 
an excitation does exist if the boundary perturbation and 
the dimensions of the cavity are favorably chosen. If the 
fractional perturbation of the boundary is given by 6 , then 
the correction terms are of that same order. 
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APPENDIX I 



• • 2 2 

It IS desirable to show that when operates 

on a function of the form cos k^x cos k^y cos k^z cos(wt+(})) , 

the result will be of the form 



2 1 

w ^ 

Q sin 



cos 



k X cos 



k^y cos k^z cos (wt+4)-3) 



where B is an angle defined such that 



1 

Q 



/ 



[1 - 

w 



tan 3 



(I.l) 



Begin by recalling that, by definition. 



c 






gt2 " Q at 



tl.2) 



Performing this operation on a function of the form 

cos k X cos k y cos k z cos(wt+(})) produces three terms, all 

with the common factors cos k x cos k y cos k z . Canceling 

X y-' z 

these like terms from both sides of (1.2) , and using the 
fact that 



k 



2 





+ k 



2 



requires that 
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2 2 ^ ^ 

-C k cos(wt+4>) + co^cos (tot+<{)) + ^ sinCwt+'J’) 

ir U 



2 

^ il^ cos(<ot+*-6) 



or, collection terms, 



OJ 



2 2 

c 2 2 , 

H - 2 ^ cos (wt+4>) + ^ sin(o)t+(f)) = ^ --r ^ - - g - cos(wt+4)-3) 



CO 



(1.3) 



Rearrangement of (I.l) gives 



[1 



c 2k2 

_E 

2 

03 



COS g 

Q sin g ' 



and substitution of this result into (1.3) requires 



2 

0 ) cos g 
Q sin g 



cos (cot+<{)) 



2 

+ ^ sin (o)t+4)) 



^ il5-e • 



2 

Canceling the common factors co /Q , and multiplying both 
sides by sin g then requires that 



cos g cos (cot +4)) + sin g sin(cot+(f)) = cos (cot+4)-g) . (1.4) 
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Using the trigonometric identity 



cosA cos B + sin A sin B = cos(A-B) 

immediately shows (1.4) to be satisfied. Thus a general 

2 2 

solution to Ql when operating on a function of a 

specific form has been shown: 



2 _ 2 



c^ PIt cos k^x cos k y cos k z cos(o)t+4)) = 

p uJ ij X y z 



2 , 

: yr COS k X COS k Y COS k„z cos (cot+(})-3) 

Q sin 6 X y'' z 



c k 

where tan 3 = p [1 - 

Q w 
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